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We study a model of dephasing (decoherence) in a two-state quantum system (qubit) coupled to 
a bath of harmonic oscillators. An exact analytic solution for the reduced dynamics of a two-state 
system in this model has been obtained previously for factorizing initial states of the combined 
system. We show that the model admits exact solutions for a large class of correlated initial states 
which are typical in the theory of quantum measurements. We derive exact expressions for the 
off-diagonal elements of the qubit density matrix, which hold for an arbitrary strength of coupling 
between the qubit and the bath. The influence of initial correlations on decoherence is considered 
for different bath spectral densities. Time behavior of the qubit entropy in the decoherence process 
is discussed. 
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I. INTRODUCTION 

There are two long-standing problems in the theory of 
open quantum systems — memory effects and the in- 
fluence of initial statistical correlations on the system 
dynamics [U |2] . These problems take on special signifi- 
cance in the theory of decoherence (i.e., the environment- 
induced destruction of quantum coherence), since the de- 
coherence time scale is usually much shorter than the 
time scales for other relaxation processes in a system. 
The well-known Nakajima-Zwanzig projection operator 
technique [31 S] and its modifications [T] provide for- 
mally exact non-Markovian master equations for almost 
arbitrary open systems and initial conditions. Unfortu- 
nately, it is impossible to solve these equations for more 
or less realistic models. Thanks to the work of many peo- 
ple (see, e.g., Refs. [H [51 [B] and references therein), we 
possess today some partial advancements in constructing 
reasonable approximations for non-Markovian dynamics 
of open quantum systems, but a systematic general ap- 
proach is still lacking. 

The simplest systems involving many of fundamen- 
tal features of quantum coherence are two-state systems. 
Such systems are important in their own right as the ele- 
mentary carriers of quantum information (qubits) [THl]. 
In addition, some two-state models admit exact solutions. 
The latter fact is very useful because it allows one to gain 
valuable insight into general properties of the dynamics 
of decoherence and can serve as a step toward consistent 
approximations for more complicated open systems. For 
instance, in Ref. [SJ [S] two exactly solvable models have 
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been used to reconstruct the corresponding exact non- 
Markovian master equations and study in detail their 
memory kernels. It would be interesting to apply a simi- 
lar approach to situations with non-negligible initial cor- 
relations between an open system and its environment. 

Recently [TD], an exact analytic solution was obtained 
for a qubit model [llH13j with initial qubit-environment 
correlations. The results demonstrate new interesting 
features of decoherence in the presence of initial corre- 
lations. Unfortunately, the initial state considered in 
Ref. [To] seems to be somewhat artificial since it relies 
on the assumption that the qubit -|- environment sys- 
tem is initially prepared in a pure quantum state at zero 
temperature. In this paper we show that the model jlll - 
IT3| admits exact solutions for a large class of physically 
reasonable correlated initial states at finite temperatures 
and derive explicit expressions for the coherences (off- 
diagonal elements of the qubit density matrix). 

The distinctive feature of the dephasing model [TIHT5] 
is that the average populations of the qubit states do not 
depend on time. In other words, there is no relaxation to 
complete equilibrium between the qubit and the environ- 
ment, i.e., the model is nonergodic. Curiously enough 
this feature may be considered as an advantage of the 
model. Indeed, since the relaxation time for dissipa- 
tive processes is usually much larger than the decoher- 
ence (dephasing) time , the model [TTHT^ corresponds 
to the limiting case Tj^/tj^ — >■ oo, when all irrelevant cor- 
rections due to energy dissipation are removed. Thus one 
may expect that essential features of decoherence in this 
simple model will be similar to those in more involved 
but less tractable dissipative models. 

The paper is structured as follows. In Sec. II we give a 
brief description of the model and find exact solutions of 
equations of motion for all relevant operators. In Sec. Ill 
we derive exact expressions for the elements of the qubit 
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density matrix, which hold for correlated initial condi- 
tions. These expressions are used in Sec. IV to study 
possible regimes of decoherence for different bath spec- 
tral densities. The results are compared with the uncor- 
related case. Finally, in Sec. V we derive an exact expres- 
sion for the qubit entropy and discuss its time behavior 
in the presence of initial qubit-environnient correlations. 



II. THE MODEL 

We consider a simple version of a spin-boson model 
describing a two-state system (S) coupled to a bath (B) 
of harmonic oscillators [1] [TTI - [T3] . In the "spin" repre- 
sentation for a qubit, the total Hamiltonian of the model 
is written as (in our units h = 1) 



with 



H = Ha 



9kbl + glhj , (1) 



where cjq is the energy difference between the excited 
state |1) and the ground state |0) of the qubit, and is 
one of the Pauli matrices ctj^ , CTj , (T3 . Note that CT3II) = 
|1) and (T3IO) = — 10). Bosonic creation and annihilation 
operators b\. and 5^, correspond to the fcth bath mode 
with frequency w^.. 

Suppose that at time t = the state of the total system 
is described by some initial density matrix g{0). Then at 
time t the density matrix is given by 

g{t) = exp {-iHt) g{0) exp (iHt) . 

Our main interest is with the reduced density matrix of 
the qubit 



0sit) = TrBe(t), 



(2) 



where Tr^ denotes the trace taken over the bath degrees 
of freedom. It is easy to see that 



{l\gs{t)\l) = ^,{l + {a,m, 
{0\gsm) = {a+it)), {l\gsitm = 



(3) 



where a^. = {a^ ± ictj) /2. Here and in the following the 
symbol {A{t)) stands for the average value of a Heisen- 
berg picture operator calculated with the initial density 
matrix of the total system: 

{A{t)) = Tr {exp {iHt)Acxp {-iHt) g{Q)} . 

The notation {A) will be used for averages dX t — Q. 

In the model ([T]), equations of motion for all relevant 
operators can be solved exactly. Technical details are 
summarized in Appendix Here we quote the results. 
The time-dependent bath operators are given by 



hit) 



(4) 



(5) 



and the qubit operators (t) can be written as 

{t) — exp [±iujQt =F R{t)] o-_|_ , (6) 
where the operator R{t) acts only on the bath states: 



R{t)=Y,['^k{t)bl~al{t)b, 

k 



(7) 



Since dg commutes with the Hamiltonian ([T]), we have 
<J^{t) = CT3, so that the populations {'^\gs{i)\^) £^nd 
(l|f?s(^)|l) do not depend on time. 



III. EXACT SOLUTIONS FOR THE 
COHERENCES 

A. Uncorrelated initial state 

As an introduction to our subsequent development, 
we first briefly review the calculation of the coherences 
((T_|_(t)) in the case that the qubit and the bath are ini- 
tially uncorrelated, and the bath is in thermal equilib- 
rium at some temperature T [3 [TT| - [T5] . In this case the 
initial density matrix of the total system is a direct prod- 
uct 



e(0) = gs{^)®gB. 



Qb 



/Zb, (8) 



where /3 = 1/k^T, and Zg is the bath partition function. 
Note that gg{0) may be a pure state as well as a mixed 
state of the qubit. 

Using expressions (|6| and ([s]), one obtains 



K(t)) = (a±)e±-o'e-^(*) 
with the decoherence function ^{t) defined as 

7(t) = - In (c^''**^)^ 



(9) 



(10) 



where the symbol {■ ■ ■) b denotes averages taken with the 
bath distribution g^. After straightforward algebra (see, 
e.g., Ref. [T]) one finds 



lit) 



duj J{lo) coth (/3w/2) ^ ^""^"^^ , (11) 



were the continuum limit of the bath modes is performed, 
and the spectral density J(w) is introduced by the rule 



duj J{uj)f{uj). 



(12) 



Expression (11) is the exact result for the decoherence 
function in the model ([!]) under the uncorrelated initial 
condition ([s]). 



3 



B. Correlated initial states 

We turn now to the coherenees {aj.{t)) for correlated 
initial states. First of all, we have to specify the initial 
density matrix g{0). Generally speaking, one may imag- 
ine a variety of different forms of this density matrix, but 
we will restrict our consideration to initial states 



that the density matrix ( 14 ) can be written as 



(17) 



where ggd'il')) plays the role of the initial density matrix 
of the bath and, at the same time, is a functional of \ip): 



g{o) = ^J2^„,e-^"nl, 



{^\exp{-f3H)\^) 



(18) 



(13) 



where operators fi„ act on the qubit states, and the par- 
tition function Z ensures the normalization of g(0) . Such 
density matrices are used, for instance, to describe the 
preparation of a system by means of a quantum mea- 
surement They also arise in a natural way in the 
calculation of correlation functions of open quantum sys- 
tems. In the theory of decoherence, of special interest 
are initial density matrices [TJ [2] 



In the model under consideration, this functional can be 
obtained explicitly. First we note that the Hamiltonian 
^ satisfies 



where 



e-^^|0) =e^"o/2e-'3<"' ® |0), 
e-^^|l) =e-^-o/2e-/3^i+' (^|1) 

Hi^^ = E ^i^^lh ± E Ubl + gib 



(19) 



(20) 



1 



P., 



(14) 



Then, with Eqs. ( 19 ), the initial bath density matrix ( 18 ) 
is manipulated to 



where P^ — is the projector onto a pure quantum 

state lijj). In the usual representation, we have 



where 



/2p-/3// 



IV) = ao|0> 



-«i|l> 



(15) 



Z^j^^ = Tr^cxp 



(21) 



(22) 



with ioqi 
written as 



= 1. 



The projector P can also be Thus, although the density matrices (|8|) and ([17| are sim 



\p\ 



(16) 



where the components of p are easily expressed in 
terms of the amplitudes Qq and ai- The density ma- 



trix (14) corresponds to the so-called selective quantum 
measurement [1] and is especially suited to study the 
environment-induced destruction of quantum coherence. 

In contrast to density matrices of the form ([S]), the den- 
sity matrices ( 13 1 and ( 14 1 contain the total Hamiltonian 



ilar in form, the latter involves initial qubit-bath corre- 
lations. As shown below, these correlations can signifi- 
cantly affect the dynamics of the total system. 



C. Time evolution of correlated initial states 

The structure of the coherences (cr_|_(t)) can be inves- 
tigated for arbitrary operators in Eq. (13 1. First we 
write 



of the system and, consequently, describe initial qubit- 
bath correlations. If we neglect the interaction term -ffj^t, 
replacing H by Hq = Hg + Hn , then we immediately re- 
cover an uncorrelated state (l8|. It is interesting to note 



ETr^Tr, [nla^n^e^"^^%-^"]} 



Using Eqs. (19), the above expression is recast into 



(-±(t)> = E {^'"°^' ml'^AlO) Tr^ (e^«We-^<-') 

+ e-^-o/2(i|j^t^^^^f^^|l)Tr^(e^H(*)e-^«r)} (23) 
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with the partition function 



It is easy to verify that 



(-) 

B 



+ e-'^-o/2(i|j^t^f^^jl)4+)| 



(24) 



We now note that the traces Tr^ (. . .) in Eq. (23) can 
be simpHfied considerably by a unitary transformation of 
i/jf ^ and R{t) with 



\9k\' 



E 



(25) 



U±R{t)U±^ = R{t) ± i$(t), 



where we have introduced the c-number function 



^ sin(w,0 



dw J{(jj) 



sinwt 



(26) 



with the same spectral density J(w) as in E q. (|TI| ). Using 
the transformation properties (251 in Eqs. (psTand (p4|, 



we obtain the final expression for the coherences, which 
can be conveniently written as 



where {(J±) are initial values of the coherences. 



Formulas ( 26 1 and ( 27|) give an exact result for the de- 
coherences in the model ([T]) with initial correlated states 
of the form (13 1. Expression ([27|) is considerably simpli- 



fied in the special case when the initial density matrix is 



given by Eq. ( 14 1. After some algebra which we omit, we 
obtain 



{a^{t)) = e^^'^otQ-iit) I cos[$(t)] 

sinh(/3a;Q/2) — (o-g) cosh(/3wo/2) . 



±i 



cosh(/?WQ/2) — (0-3) sinh(/3wQ/2) 



sin[$(i)] ,(28) 



where (ug) = joj^p — jagp. This expression can be repre- 
sented in a more transparent form 

(cr±(t)) = ((T±) e±*["o*+x(t)le-7(t)^ (29) 



where 



7W=7W+7corrW 



(30) 



is the modified decoherence function which includes the 
correlation contribution 



7corr(0 



In 



[cosh(/3wo/2) - (ag) sinh(/3wo/2)] 
and x(t) is the time-dependent phase shift with 
sinh(/3ajo/2) — (^3) cosh(/3wo/2 



,(31) 



cosh(/3a;Q/2) — (cTg) sinh(^Ci;g/2) 



tan[x(i)] = l^v,;^. /o^ ,. \Zl^)o tan[$(t)]. 

(32) 

It is a straightforward matter to derive a formula analo- 
gous to Eq. ( 29 1 for more general initial states ( 13 1 be- 



cause each of the operators in Eq. (271 can always 



be represented as = + a • with some (in gen- 
eral, complex) and p^. The resulting expressions for 
7corr(0 x{i) a-i'G rather cumbersome and will not be 
given here. From now on we restrict our discussion to 
Eqs. ([29|-([32l). 



IV. REGIMES OF DECOHERENCE FOR 
DIFFERENT BATH SPECTRAL DENSITIES 

As already noted, the model ([I]) is noncrgodic, i.e., 
it does not describe the establishment of complete ther- 
mal equilibrium between the qubit and the bath since 
((Tg(i)) — const. In this connection it is of interest to in- 
vestigate the long-time behavior of the decoherence func- 
tion (30). Suppose, for instance, that '^{t) is a bounded 
function. Then, in the limit as t — > cx), the averages 
{(j^{t)) do not tend to zero, so that the destruction of 
quantum coherence is "incomplete" . Physically, in this 
case the final stage of decoherence is determined by slow 
processes involving exchange of energy between the qubit 
and the environment, which are not included into the 
Hamiltonian ([ll). On the other hand, if ^{t) — ^ cxi as 



t — > 00, then (o\j_(f)) — > showing complete decoherence 
within the framework of the model under consideration. 



Note that the correlation corrections in Eq. (31 1 [or 



for a more general case, in Eq. (27)] are always bounded 



This means that initial correlations alone cannot lead to 
complete decoherence. Let us now turn to the "dynami- 
cal" part of the decoherence function given by Eq. (111. 
It can be written conveniently as a sum 



7(0 =7vac(i) +7th(*)> 



(33) 
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where 



7vac(*) 



duj J{lo] 



1 — cos uit 



(34) 



is the contribution to the decoherence function from vac- 
uum fluctuations in the bath, and 



7th (i) 



dujj{uj) [coth (/3w/2) - 1] 



1 — cos Ujt 



duj 



J{uj) 1 — coscji 
- 1 Z/^ 



(35) 



is the temperature dependent thermal contribution. We 
have obvious inequahties 



7vac(i) < 2 

7th (i) < 4 



duj 





oo 



J{u) 



w2 (e/3<^ _ 1) 



(36a) 



(36b) 



which show that 7vac(0 ^-^"^ 7th (^) ^'"'^ bounded if the 
integrals converge. 

Another important question is whether, and under 
what conditions, the total decoherence function 7(t) 
tends to a definite limit as i — > cxd. Necessary conditions 
for different terms in Eq. ( 30 1 to have definite long-time 



limits can be derived by using relation |14j 



lim f{t) = lim £ / dte-'*f{t), (37) 

t-i-oo Jq 

which is valid if the limit on the left-hand side exists. 



When apphed to functions ([Sff, (|35|, and (126|), Eq. (jST]) 

gives 



lim 7^ 

t—^oo 



At) 



du 



J{UJ) 



nOO 

Hm 7tjj(i) = 2 / duj 



J(w) 



,2 _ 1) 



lim $(t) = - lim -^-^ 

i— >oo 2 >0 io 



(38a) 



(38b) 



(38c) 



To go beyond relations ( 36 ) and ( |38[ ) , one needs some 
information about the bath spectral density. Generally 
speaking, its form can be obtained from a fully macro- 
scopic analysis of the system-bath interactions leading 
to the spin-boson model ([l]). It would be beyond the 
scope of this paper to present a detailed discussion of 
such derivations. In many cases of practical interest (see, 
e.g., Refs. |13|[15]), J(w) may be considered to be a rea- 
sonably smooth function which has a power-law behav- 
ior J(w) oc w** (s > 0) at frequencies much less than 
some "cutoff" frequency fi, characteristic of the bath 
modes. In the limit w — >■ oo, J(w) is assumed to fall 
off at least as some negative power of lo. Then, on di- 
mensional grounds, we may write 



j{uj) = x^n^-'^Lu'' F{Lu/n), 



(39) 



where is a dimensionless coupling constant, and a cut- 
off function satisfies 



F(0) = 1, 



lim 



F{uj/n) = 0. 



(40) 



The case s = 1 is usually called the "ohmic" case, the case 
s > 1 "superohmic" , and the case < s < 1 "subohmic" . 

As is easy to see, the convergence of the integrals in 
Eqs. ( 36 1 and ( 38 1 depends crucially on the low-frequency 
behavior of J(w). Assuming a power law for J{u:) in 
this region, we can make some conclusions about possible 
regimes of decoherence in the model ([I]). Let us first of 
all note that both integrals in inequalities ( 36 1 converge 
for s > 2. Therefore, in this case the total decoherence 
function ^{t) is bounded, so that the above mentioned 
regime of "incomplete decoherence" takes place. It is 
also seen from Eqs. ( [36| that, in the superohmic case 
with 1 < s < 2, the thermal term 7tij(i) is the only 
contribution to the total decoherence function which may 
diverge as t — oo. Finally, Eq. (38c) shows that the 
correlation term 7corr(^) may have a nonzero long-time 
limit only in the ohmic case (s = 1). 

The virtue of the above conclusions is that they apply 
to all spectral densities of the form (39), but they say 



nothing about the regime of decoherence for < s < 2. 
It is reasonable to expect that in this range of values of 
the parameter s we are dealing with the regime of "com- 
plete decoherence" , since the vacuum term 7^^^ [t) may 
diverge as i — !■ 00 for < s < 1 and the thermal term 
7th (0 for < s < 2. We could not determine, how- 
ever, a general sufficient condition for the cutoff function 
in Eq. (39) to ensure that 7vac(*) ^"^^ 7th (^) behave in 



this manner. Note, however, that this question is largely 
academic since the model specified by the Hamiltonian 
([1]), on its own, is of most physical interest in the low- 
temperature range (i7/3 ^ 1) where dissipative processes 
are irrelevant and the bath dynamics is completely de- 
termined by low- frequency modes with Wj, ^ f7. For this 
temperature range, we show in Appendix [B] that, under 
assumptions (39) and (40 1, the thermal term 7th(i) is a 
steadily increasing function of time for < s < 2, and 
hence 7th (^) — >■ 00 as t — >■ 00 . This means that the regime 
of "complete decoherence" takes place for all < s < 2. 

Up to this point we have been concerned with those 
features of decoherence which do not depend on the form 
of the bath spectral density or, in particular, on the form 
of the cutoff function in Eq. ([39|). For specific cases, func- 



tions (31), (34), and (35) can be evaluated numerically. 



However, to gain more insight into the overall picture of 
decoherence, it will be instructive to consider a particu- 
lar choice of the cutoff function in Eq. ( 39 ) for which all 
the quantities of interest can be evaluated exactly. We 
shall take for the bath spectral density the expression 



-uj/n 



(41) 



which is most commonly used in the theory of spin-boson 
systems [HinHIllIi]. 
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Let us start with the vacuum term (34) in the deco- 



herence function. Substituting here the spectral density 
from Eq. (|4T]) and doing standard integrals, we get 



7vacW = A,r(s-l) 

J cos [(s - 1) tan-i(m)] 1 

""l' (l + f}2i2)(-l)/2 J' 



(s^l), (42a) 



(42b) 



The latter expression is identical with the well known 
result p]. Note that the properties of 7vac(*) given by 
Eqs. (42) are consistent with the general statements for- 



mulated in the beginning of this section. First, the vac- 
uum term is a monotonically increasing function of time 



for s < 1. Second, in the superohmic case (s > 1) this 
term has a long-time limit: 



lim7,,,(i) = A,r(s-l), (s>l). 



(43) 



It is easy to check that this limit is exactly the same as 
given by Eq. ( 38a ) . One can also see from Eq. ( 42a I that 
7yjjj.(i) monotonically saturates to j^^^{oo) for 1 < s < 2 
and is a nonmonotonic function of time for s > 2. These 
properties of the vacuum term are illustrated in Fig. [T] 



Evaluation of the thermal term ( 35 ) with the spectral 



density ( 41 ) also reduces to performing standard inte- 



grals. After some manipulations which we omit, we ob- 
tain 



7th(0 = 2A, m'-' Vis - 1) 



OO 

V - 



{k + i/nfiY 



-(s-l)/2 



7th (0 = 2Ai 



InF (1 + l/np) - ^ In |r (1 + l/np + it/l3)f 



cos[(s-l)<^fe(t)]|, (s>0, s^l), (44a) 
(s = 1), (44b) 



where we have introduced the notation 

.-1 / t/P 



fkit) = tan" 



fc + 1/0/3 



(45) 



where 



Tg — /S/tt = h/Trk^T 



(47) 




is the so-called thermal correlation time. 



FIG. 1. Time dependence of the vacuum contribution to the 
decoherence function in the superohmic case (s > 1). 



At low temperatures (fi/? ^ 1) expression (44b) re- 
duces to the well known result pQ O [T^] 

7th(<)--Ailn|F(l + »i//3)|' 
sinh(t/r^ 



Ai In 



(s = 1, l]/3> 1), (46) 



We see from Eq. (44a) that in the case of incomplete 



decoherence (s > 2) the thermal term has a long-time 
limit 



lim 7,h(i) - 2A, m'-' 

I— J-OD 

X T{s -l)C{s- 1,1 + l/ni3), (s>2), (48) 



where ^{z,v) is the generalized Riemann zeta function. 
One may easily verify that the above result is exactly the 
same as predicted by Eq. ( 38b I . It is somewhat more dif- 



ficult to clarify the time behavior of 7th(i) directly from 
Eq. (|44|). We show in Appendix Ib| that, for < s < 2, 



7^jj(t) is a monotonically increasing function of time and 
7th (^) ^ OO as i OO. In the regime of "incomplete 
decoherence" (s > 2), the thermal term exhibits non- 
monotonic time behavior for sufhciently large values of 
the parameter s (see Appendix [b| and Fig. [2]). 



We now turn to the correlation term (31) in the deco- 



herence function assuming the bath spectral density to 
be given by Eq. (41 1. As a preliminary step, we calculate 
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FIG. 2. Time dependence of the thermal contribution to the 
decoherence function in the regime of "incomplete decoher- 
ence" (s > 2); Tg is given by Eq. (47 1 and firg — 10. 



the phase function ( |26[ ). Doing the cj- integral, we find 
\r(s-l) 

- (l + f^2,2)(-l)/2 

X sin [{s - 1) tan~i(m)] , (s > 0, s ^ 1), (49a) 

= Ai tan-i(m), (s = 1). (49b) 

These expressions allow one to bring out some important 
properties of the correlation contribution to the decoher- 
ence function. First, in the subohmic case (0 < s < 1), 
the phase function (49a) increases with time and, conse- 



quently, 7corr(0 oscillates. Second, for s > 1 the phase 
function has long-time limits 



lim $(t) = Ai7r/2, (s = 1), 



lim $(<) 0, 



(s>l). 



(50) 



which are consistent with Eq. (38c I. We thus conclude 



that in the superohmic case (s > 1) the correlation term 
(311 asymptotically tends to zero as t — oo, whereas in 
the ohmic case (s = 1) it has a long-time limit 



lim 7corr(0 



(1 



sin^(Ai7r/2) 



[cosh(/3wo/2) - (0-3) sinh(^Wo/2)]' 



(51) 



It is interesting to note that this limiting value is a peri- 
odic function of the coupling constant Aj^. In particular, 
for Xi — 2n, (n — 1,2,.. .), we have 7(.o„(oo) — 0. The 
time behavior of 7corr(0 i'^ the ohmic case is very sensi- 
tive to the value of the coupling constant. For sufhciently 
large A^, 7coit(0 h^-s a "peak" structure (see Fig. [3]). 

Qualitative properties of the contributions to the total 
decoherence function 7(f) for different values of the pa- 
rameter s in Eq. (41 ) are summarized in Table|l] We note 



FIG. 3. Time dependence of the correlation contribution to 
the decoherence function in the ohmic case (s = 1) for dif- 
ferent values of the coupling constant A^. Parameter values; 
(as) = (equal populations of the qubit states |0) and |1)); 
a;o/3 = 1. 



that the value s = 2 plays a role of a "critical parame- 
ter" for decoherence in the model ([ij. If s < 2, we have 
the regime of "complete decoherence" since 7(t) — > 00 
as t — > 00, and hence the coherences ([29]) asymptotically 
tend to zero. For s > 2, the total decoherence function 
has a finite long-time limit determined by the vacuum 
and thermal contributions. In this case we are dealing 
with the regime of "incomplete decoherence" . It should, 
however, be pointed out that in a real qubit the resid- 
ual coherences (a_|_(oo)) decay to zero due to dissipative 
processes which are not included in the model ([T]). 

We close this section with remarks about the role of 
initial qubit-bath correlations in different regimes of de- 
coherence. Some conclusions concerning correlation ef- 
fects can be drawn directly from Eq. (31). First, in the 
weak coupling limit (A^ ^ 1) Eqs. (49) and ( [3l] ) give 
7j.o;.;. oc A^, while ^^^^ cx A,, and °^ '^s- We see that in 
this limit the main contribution to the total decoherence 
function ( [30| is from its "dynamical" part 7(<). Second, 
it is clear that the correlation term (31) is small com- 



pared to 7(t) at extremely low temperatures (/Swg 3> 1) 
for all values of A^,. We thus conclude that the role of ini- 
tial qubit-bath correlations becomes pronounced in the 
temperature region Ploq < 1 for intermediate strength of 
coupling. 

Formulas (49) show that $7"^ determines the charac- 
teristic time scale for the correlation effects. The same is 
true for the vacuum term (42). On the other hand, the 



quantity given by Eq. (47) determines the time scale 
for thermal effects in decoherence [cf. Eqs. (44) and (45)]. 
In the most interesting low-temperature range (J7/3 3> 1), 
we have fir^ ^ 1, so that initial correlations might be 
expected to have a pronounced effect on the coherences 
(cr_|_(i)) at times t < Tg when the main contribution to 
the "dynamical part" ^{t) of the decoherence function 
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TABLE I. Time behavior of different contributions to the total decoherence function "/{t) for the bath spectral density (41 1. 



7th (t) 



7c orr it) 



< s < 1 
s = 1 

1 < s < 2 
s > 2 



Monotonic increase 
Monotonic increase 

Saturation to 7vac(oo) / 
Saturation to 7vac(oo) 7^ 



Monotonic increase 
Monotonic increase 

Monotonic increase 
Saturation to 7^1^(00) 7^ 



Oscillations 
"Peak" structure; 

7corr(00) = Or 7corr(oo) / 

depending on the value of 
Nonmonotonic decay 
Nonmonotonic decay 



comes from vacuum fluctuations in the bath. For the su- 
perohmic case, this can be clearly seen in Fig. |4] which 
illustrates the time behavior of the coherences in differ- 
ent regimes. Note, however, that in the subohmic and 
ohmic cases, for sufficiently weak coupling, the correla- 
tion effects manifest themselves at times t > Tg where 
thermal excitations dominate (see Fig. |4]). The reason 
can be traced to the facts that, in the former case, the 
correlation term 7coit(^) exhibits undamped oscillations, 
while in the latter case it has a long-time "plateau" . 



V. ENTROPY 

Entropy plays a crucial role in the theory of open sys- 
tems since it is a natural measure of the lack of informa- 
tion about a system. It is thus of interest to discuss the 
time behavior of the qubit entropy for the model under 
consideration. 

We start from the general expression for the von Neu- 
mann (information) entropy of a quantum system: 



S{t)^~Tvs{gsit)\ngs{t)}. 



(52) 



In the case of a qubit, it is convenient to express S{t) in 
terms of the Bloch vector v{t) — {(j{t)) using two rep- 
resentations for the density matrix gg{t). The standard 
representation is [TH [T7| 



(53) 



In Appendix [C] we derive another representation which 
is better suited to calculate \nQg{t): 



(t) - I ^l~v^t) exp [a ■ uit)] , (54) 



where u X] v, and 



In 



1 



V 



(55) 



Strictly speaking, formula (54) is valid only if w < 1, 
i.e., for a mixed state. Note, however, that the limit 
V —?' 1 can be taken directly in the entropy (52) after 



calculating the trace. With expressions (54) and (55) 
one easily derives from Eq. (52) 

1 



(l-w)ln(l-t;) 



(56) 

For a pure state {v — > 1) this formula gives 5 = 0, as it 
should be. 

The square modulus of the Bloch vector can in general 
be written as 

v\t)=4{a4t)){a_it)) + {a,{t)f. (57) 

For simplicity we shall restrict further discussion to cor- 



related initial states of the form (14). Then we have 
v{0) = 1 and, consequently. 



Now using the solution (|29|) and taking into account that 

(58) 



CT3 is an integral of motion, we obtain from Eq. ( 57 ) 

1/2 



-27(t) 



Formulas ( 56 ) and ( 58 ) determine the time evolution of 



the qubit entropy. 

In discussing the properties of entropy in the model ^ 
under the assumption ( 39 ) for the bath spectral density. 



it is necessary to distinguish two cases: the regime of 
"complete decoherence" , and the regime of "incomplete 
decoherence" . In the former case (s < 2) we have j{t) — 
00 as i — cx), and hence the limiting value of the entropy 
is the same for both (correlated and uncorrelated) initial 
conditions: 

I(l + |(a3)|)ln(l + |(a3)|) 



lim S{t) 

t^OO 



In 2 



3)|)ln(l-|(a3)|), (.<2). (59) 



The maximum entropy S^^g^^{oo) = ln2 corresponds to 
the initial state with equal populations {{cr^) =0). A 
similar situation occurs in the case of "incomplete deco- 
herence" (s > 2) since 7corr(*) — >■ as i — >■ 00. Note, 
however, that the limiting value of the qubit entropy is 
now given by 

1 
2 



lim S{t) 

t—>-oo 



In 2 



(1 



1 



ln(l + «co) 
), (s>2), 



(60) 
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1 - y, A=0.2 

2-y. A=0.2 

3 - y, A=l 

4-y, ,1=1 



1.0 

0.8 
0.6 
0.4 
0.2 
0.0 



Kcr«>| 



1 - y, A=0.5 
2-y, A=0.5 
3 - y, A=3 
4-y, A=3 





0.0 0.2 0.4 0.6 



FIG. 4. Time dependence of the coherences {^{t)) = {a±{t)) for the correlated initial condition (141 (solid lines) and the 
uncorrelated initial condition (|8| (dashed lines). In all cases A = A^; is given by Eq. (471. Parameter values: O.Tb = 10, 
cjQ^S = 0.1. 



where 



i'^sf + (1 - Ws)') 



-27(00) 



1/2 



(61) 



Recalling Eqs. (43) and (48 1, we have 



7(c^) = A,r(s-l) 
2 



1 



jC{s-i,i + i/m 



(s > 2). (62) 



The second term in braces corresponds to the contribu- 
tion from thermal excitations in the bath and is relatively 
small in the low-temperature range {ft(3 3> 1). 

Although initial qubit-bath correlations do not con- 
tribute to S{oo) for all s > 0, they influence the behavior 
of the qubit entropy at times t < Tg. The main reason is 
easy to see when one recalls that, for sufficiently strong 
coupling, the modulus of the Bloch vector v{t) may have 
rather sharp peaks associated with the evolution of initial 
correlations (see, e.g.. Fig. |4]). 

Especially interesting is the ohmic case (s = 1) where 
the time dependence of the coherences is very sensitive 
to the value of the coupling constant A^. Figure |5] illus- 
trates the kind of the behavior of the qubit entropy one 
might expect in this case. We call attention to the fact 




0.00 



0.05 



0.10 



0.15 



0.20 



FIG. 5. Time evolution of the qubit entropy in the ohmic case 
for different coupling strengths. Parameter values: firg — 10, 
LJol3 = 0.1. 



that in some time intervals the entropy for a larger cou- 
pling constant is lower than that for a smaller coupling 
constant. 



VI. CONCLUSIONS 

Let us present a summary of the results obtained in 
this work and discuss their relation to problems in open 
system dynamics. 



We have derived exact formulas (27 1 and (281 which 
describe the time behavior of the coherences (a^t)) (off- 
diagonal elements of the qubit density matrix) for all 
correlated initial states of the form (13) and ([T4|, re- 



spectively. Initial states of these types can be interpreted 
physically in terms of the general theory of quantum mea- 
surements [H IT^I^U] . In particular, operators enter- 
ing Eq. (27) represent operations during the "prepara- 
tion" of a combined qubit-bath system in an initial state 
at temperature T = 1//3. 

In this paper we have studied in detail the reduced 



qubit dynamics for the initial condition ( 14 ) which cor- 



responds to the preparation of the qubit in some pure 
state Itp). Since the qubit is not isolated from the bath, 
initial correlations come into play through the initial bath 
density matrix (18) which is a functional of \tp) [see also 
Eq. 



The quantity of most interest is the decoherence func- 
tion j{t) which determines the decay of the coherences: 
(cr_|_(t)) oc exp[— 7(t)]. Within the framework of the 
model under consideration, the decoherence function can 
be split into the vacuum part 7vac(0i the thermal part 
Jth{t), and the "correlation" part 7corr(^)- Sec. IV we 
have established some connections (with different levels 
of generality) between the properties of the bath spec- 
tral density J{u)) and possible regimes of decoherence. 
First, we have simple but quite general relations (36) 



and ( 38 ) which show that the long-time behavior of all 



the contributions to the total decoherence function de- 
pends crucially on the form of J(w) in the low-frequency 
range. In particular, inequalities (36) give sufficient con- 



ditions for the regime of "incomplete decoherence" when 
{cFj.{t)) do not tend to zero as t — > 00. Second, assum- 
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ing only that in the low- frequency range J{uj) scales as 
J(a;) oc with s > 0, we have shown that the regime of 
"incomplete decoherence" occurs for s > 2 when the to- 
tal decoherence function is bounded. Finally, taking the 
bath spectral density in the form ( |41[ ), we have derived 
exact expressions for all the terms in 'y{t). As expected, 
these expressions confirm the earlier general predictions. 

We have seen that the qualitative features of decoher- 
ence are essentially determined by the behavior of J{uj) 
in the range a; ^ 57, where is a cutoff frequency, char- 
acteristic of the bath modes. Thus, although explicit re- 
sults have been obtained for a special (exponential) form 
of the cutoff function in J{uj), there is a good reason to 
think that the conclusions summarized in Table |T] apply 
as well to other forms of the cutoff function. 

Decoherence is a fundamental property of open quan- 
tum systems which are met in quantum optics, electron- 
ics, atomic and molecular physics, etc. The literature 
concerning diverse aspects of this phenomenon is now 
quite voluminous. At first sight a study of very simple 
models like that specified by the Hamiltonian ([T]) might 
appear to be essentially a problem of applied mathemat- 
ics. In a sense this is true. It should be noted, however, 
that the pure dephasing mechanism of decoherence de- 
scribed by the model ([T]) can dominate in real physical 
systems [U] [55] . From this point of view exactly solvable 
dephasing models seem as themselves to deserve thor- 
ough studies. Note also that exact results obtained for a 
dephasing model can serve for constructing approximate 
solutions for more complicated dissipative systems [T| if 
the dissipative coupling with environment is weak. An- 
other field of applications of exactly solvable dephasing 
models with initial system-environment correlations is 
the theory of quantum information. As a first step it 
would be desirable to generalize the results to many-qubit 
systems (quantum registers). 

Recent works [101 US] show that initial system- 
environment correlations can play an interesting and 
somewhat unexpected role in decoherence. We hope that 
the results of the present paper will be useful for further 



studies of correlations effects in decoherence phenomena. 
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Appendix A: Heisenberg equations of motion 

Since CTg commutes with the Hamiltonian ([T]), equa- 
tions of motion for 6^. {t) and 6^, (t) are 

dbJt) , , , dbUt) , 

The solution of these equations with initial conditions 
&j,(0) = hf. and b\(Q) — b\ is easily found to be given by 
Eqs. Q. Equation of motion for cr_|_(t) reads: 



da,{t) 

i ; 

dt 



3a+(t) - 2^ [g^blit) + glb^it)] a^t). 



Using Eqs. Q and relation a^a_^_{t) — <T_^_{t), we get 
.d(7,{t) 



dt 



(Al) 



with 



W{t)^2Y,{9ke"^^'[bl+al{t)/2] 

k 

+ 5^6— '=*[5fe-fa,(i)/2]}. 



The solution of Eq. ( Al ) is given by 



cr_i_(t) — exp (iujQt) exp 



+ 



i / W{T)dT 







' + ' 



(A2) 



(A3) 



where exp_|_ [. . .] is the chronologically ordered exponent. 
Expression ( A3 1 can be simplified using the operator 
identity 



exp. 



I I A{t) dr 





exp 



dri / dT^ [A{t^),A{t^)] \ exp 



I I A{T)dT 




(A4) 



which is valid if the commutator [yl(rj^), ^(Tg)] at differ- 
ent times is a c-number function. It is easily verified that 



the operator ( A2 ) satisfies 



Applying Eq. ( A4 ) to W and evaluating the integrals on 
the r.h.s. of that formula, one shows that the c-number 
functions cancel. Then from Eq. ( A3 ) follow expressions 
S for (j^[t) and a_{t) = a\{t). 



[W{t,), W{^)] = -8*^ sinK(ti - t^)]. 

k 
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Appendix B: Time derivative of 7th (*) 

Let us consider the time derivative of 7tii(i) given by 
Eq. (pSl). We have 



dt 



= 2 



duj 



J{uj) sinujt 



_ 1 



Under the assumption ( 39 1 , this reads 



(Bl) 



where <fj.{t) is given by Eq. (45). Since < (pk{t) < n/2 
for all t > 0, we may employ the same arguments as 
above. Thus, for < s < 2, the thermal term (0 
monotonically increases with time and 7 ^^ (t) — > cx) as 
t — > 00. A more elaborate analysis of Eqs. ( |B3[ ) and ( |B5[ ), 
which will not be given here, shows that in the regime of 
"incomplete decoherence" the derivative d'y^y^{t)/dt is a 
positive definite function of time for 2 < s < 3 and may 
change sign for s > 3. In the latter case the thermal term 
7tjj(i) exhibits nonmonotonic time behavior. 



d7th 2A, (17/3)1- 



dt 



dx 



sm{xt/f3) . 



(B2) 

In the limit n/3 3> 1, the cutoff function may be replaced 
by F{0) = 1. Then, performing the x- integral, we get 



dj,^ ^ 2A,(»/3)i 
dt P 



E 

k=l 



n [s tan^i(t//3fc)] 



21^2 



(B3) 



where r(s) is the gamma function. If < s < 2, we see 
that dj^Y^/dt > for any t > 0. However, this fact is 
not sufficient to conclude that the decoherence function 
diverges as f — )■ cx) since it does not exclude the possi- 
bility that d'jth/dt — ?> 0+ as t — » 00 and 7tjj(i) saturates 
to a finite value. Relation (38b) excludes this possibility. 
It tells us that, for < s < 2, the decoherence func- 
tion "f^Y^it) cannot have a definite long-time limit since, 
if it had, the integral on the r.h.s. would be finite. We 
thus conclude that 7th (0 ^ monotonically increasing 
function and tends to infinity as t ^ 00. 

If the bath spectral density is given by Eq. ( [4l] ), the 
cutoff function in Eq. (B2) has the form F{x]Tlj3) — 
exp {—x/Vlf5). Now the integration over x leads to 



d7th(0 _ 2A, (17/3)1 



dt 



■ ns) 



sin [s if kit)] 



fc=i [(fc + i/npf + [t/pY 



s/2 ■ 



(B4) 

{s > 0), (B5) 



Appendix C: Exponential form of the density matrix 



The derivation of formula ( 54 ) is based on the following 



property of the Pauli matrices (see Problem 1 to § 55 in 
Ref. [E]). Let / be some function of the operator a+a-u, 
where a and u are c- number quantities. Then 



where 



f{a + a ■ u) = A + (f ■ B, 



B 



1 u 

2 u 



[f{a + u)-f{a^u)]^ 



(CI) 



(C2) 



and u = \u\. Identifying the r.h.s. of Eq. (CI) with the 
qubit density matrix written in the form ( 53[), we have 
A = 1/2 and B = v/2. On the l.h.s. of Eq.^H| we take 
f{x) = exp{x). Then Eqs. (IC2]) give 



1 = e" (e" + e"") , u = e'^ (e" - e" 
Solving these equations for a and u, wc find 



1 1 



l-v 



This leads immediately to the representation ( 54 ) for the 
density matrix. 
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